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We analyze the influence of electron-positron pairs creation on the motion of vortex 
lines in electromagnetic field. In our approach the electric and magnetic fields satisfy 
nonlinear equations derived from the Euler-Heisenberg effective Lagrangian. We 
show that these nonlinearities may change the evolution of vortices. 
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I. INTRODUCTION 

The phenomena of creation and evolution of vortices have always attracted people's at- 
tention both in the past and in the present. In contemporary physics they gained particular 




interest since having been experimentally observed in Bose-Einstein condensates [1|,|2, 
Vortices in superfiuids, due to the absence of viscosity, exhibit certain unconventional fea- 
tures like the persistence of the whirl or its singular nature. The Bose-Einstein conden- 
sate may be described by the nonlinear Schrodinger equation (the Gross-Pitaevskii equa- 
tion satisfied by certain macroscopic wave function. In that way one has been led 
to studying vortices in quantum mechanics (QM). There is a striking resemblance between 
dynamics of fluids and QM via the hydrodynamic formulation of the latter [8] . QM can, even 
in the linear version, serve as a model theory for investigating the behaviour of vortices in 
superfiuids. Such studies, concerned with the dynamical as well as 



tex evolution in various configurations, both in nonlinear 



and in linear 
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.s well as topological aspects of vor- 

m 



23| cases have recently been undertaken. 



Together with the attention paid to nonrelativistic QM the singular solutions in other 
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field theories as electromagnetism for instance liave been investigated 
will be concerned with this question also in the present paper. While considering vortices 
in fields corresponding to spinning particles one encounters the problem, that the wave 
funcion has more than one component and the condition ip{r^t) = 0, leads to too many 
equations which cannot be simultaneously satisfied. One of the solutions of this problem for 
electromagnetic field was proposed in Ref. |27|, where vortex lines were defined by the null 
values of two relativistic invariants: S and V. These two equations mean two surfaces in 
the three-dimensional space. Their intersection in general may be a curve — a vortex line. 
We will base our work on this approach. 

In the present work we would like to investigate how quantum effects can influence the 
motion of the nodal lines of the electromagnetic wave function [2^ • The term "quantum" is 
used here in the field theoretical sense: Maxwell electrodynamics, as well as Schrodinger wave 
mechanics, are classical from that point of view. We will concentrate on the effect of electron- 
positron pairs creation in electromagnetic fields. In order to construct the position and time 
dependent wave function, we still need the classical equations for electric and magnetic 
fields. These equations are no longer linear since pairs creation leads to the photon-photon 
interaction and Maxwell equations are suplemented by additional terms which, in the lowest 
approximation, are cubic in fields and quadratic in the fine- structure constant (in general 
they might be also nonlocal). Although the correction is small it can certainly influence the 
motion of the vortex lines and particularly change their topology. 

One should mention here that there exists also another type of quantum effects — which 
remain beyond the concern of the present work — connected not with the e~^e~ content of 
the vacuum, but with fluctuations of the electromagnetic field itself. These effects, due to 
the nonzero vacuum expectation value of billinears in fields, lead to the smoothing of the 
vortex core (i.e. the line on which both invariants are equal to zero). In this case the core is 
no more singular. It is defined not by the condition F"^ = 0, which is not satisfied, but rather 

^ , where ~ means |-^|c/assica« < l^llac.fiuct. S]; where F is the Riemann-Silberstein 
vector spoken more of in Sections II and III. 
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describing. 



As a starting point we choose the Euler-Heisenberg (E-H) Lagrangian 
in the lowest order, the dynamics of classical electromagnetic fields with vacuum polarization 
effects taken into account. Field equations obtained from this effective theory in Section II, 
exhibit solutions containing vortex lines, the evolution of which may be viewed and compared 
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to that obtained from the classical Maxwell equations. In this work we analyze two such 
cases. Both are chosen from to make the comparison of the results in our works very 
easy. Our results are presented in Section III. 

The main practical problem in this investigation comes from the fact that quantum 
corrections are, in general, small and it is very hard to see them on a drawing. The choice of 
examples considered in our work from among those of Ref.^ is dictated just by the criterion 
of quantum effects being visible. It is clear that they are noticeable not by analyzing or 
measuring the precise shape of a vortex line — the slight deviation of which from that 
obtained in classical theory surely does occur — but rather by observing "to be or not to 
be" effects or topological ones. 

There are two limitations which cause some of the results of this work to be qualitative 
rather than quantitative. Firstly they are obtained within perturbative regime. This regime 
means that the electromagnetic field may not be too strong and its strength is limited by 
the condition of a\F\'^/m^ being small. Fields considered here are polynomial, so this re- 
quirement means that the evolution should not go beyond certain limited space-time region. 
However, close to the vortex line, similarly to the situation in Gross-Pitaevskii equation, the 
terms quadratic in S and V in the Hamiltonian © become small, even for large values of 
electromagnetic fields, and perturbative calculation is again well justified. 

Secondly we have to remember that that E-H Lagrangian describes only slowly varying 
fields, for which the nonlocality may be neglected. Their relative change at a distance of 
the Compton wavelength of the electron should be small. In view of that the E-H effective 
Lagrangian is treated in our work as a certain nonlinear model of the true theory of elec- 
tromagnetic fields obtained form QED without real charges. An another interesting model 
in this context constitutes the Born-Infeld electrodynamics js^. One should, however, have 
in mind that even small corrections, coming from weak fields can change the topology of 
vortices. 

II. FIELD EQUATIONS 

The Euler-Heisenberg Lagrangian Isil. Is^. which accounts for the vacuum polarization 
processes in the lowest approximation, has the following form 

C{r, t) = S{r, t) + ^ [4S{r, tf + 7P(r, tf] , (1) 
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where S and V denote the two Poincare invariants formed of electromagnetic fields 

S=^{E^ -B^) , V = E B (2) 

and a and m are fine-structure constant and electron mass recpectively. The canonical 
variables are electric and magnetic inductions: D and B, where the former plays the role 
of canonical momentum and the latter of position In this picture the electric field 

strength E in the Lagrangian corresponds to velocity. The field equations, we will need 
for our purpose, are the canonical Hamilton equations 

where 7i(r,t) denotes the Hamiltonian density. To find the explicit form of Q we have to 
perform the Legendre transform and pass from C to Ti. The canonical momentum is, as 
always, defined as a derivative of the Lagrangian over velocity 

Dir t) = ^^("'^^ = ^-^E(r t) + t) (A) 



1 + jSir^t) 



E{r^) + ^,r{r,t)B{r^). (5) 



which gives 

This kind of equation usually bears the name of a constitutive equation. It reflects the 
nontrivial structure of the medium. In the present case this medium is the quantum field 
theory vacuum with its polarizability via electron-positron pairs creation and anihilation. 
We now need to invert this equation and express velocity E in terms of canonical variables 
D and B. Since our initial Lagrangian is given only in one loop approximation (a^) 
then our further calculations may be led up to this order too. We can therefore postulate 
E{r,t) in the form 

E{r, t)=[l + a^}C{r, t)] D{r, t) + a'Mir, t)B{r, t) , (6) 

where quantities /C(r,t) and Ai{r,t) are to be determined. Substituting (jH)) into ne- 
glecting terms of the order higher than a^, and comparing coefficients multiplying vectors 
D and B we find 

1 fi 

/C(r,t) = -^[i?(r,tf-B(r,tf] , (7a) 
Mir,t) = -^D(r,t).B(r,t). (7b) 
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The Hamiltonian density may be now found as 



n{r,t) = E{r,t) ■ D{r,t) - C{r,t) , 



(8) 



where E in the whole above expression should be eliminated in favor of D and B, according 
to the relations ()6l7a|) . The explicit form of Ti. is then 



1 



2a' 



[D{r,tf~B{r,t)'] 



212 



14^2 



2 L ^' ■ ^' J 45m4 L ^ ' ^ ^ ' ^ J 45m4 
Now we are in a position to write down the equations (jH)) in an explicit form 



[D{r,t)-Bir,t)f . 

(9) 



D{r,t) = V X <^ B{r,t) 



{D{r,tf-B{r,tf) 



(10a) 



45m4 



D{r,t) [D{r,t)-B{r,t)] 



B{r,t) = - V X <^ D{r,t) 



1 



45m4 
28^2 



{D{r,tf~B{r,tr) 



(10b) 



45m^ 



S(r,t) [£)(r,t)-S(r,t)] 



Introducing two complex vectors F±{r,t) according to the relation 

F^{r,t) = ^{D{r,t)±tB{r,t)) 



(11) 



we can rewrite the equations (|Tn|l in the form 

F4r,t) = -zVxF+(r,t) + ^Vx [F_(r, t) (llF+(r, t)^ - 3F_(r, t)^)] , (12a) 
F_{r,t) = ^VxF_(r,t) + |^Vx [F+(r, t) (llF_(r, t)^ - 3F+(r, t)^)] . (12b) 

In the classical case the right hand sides of ()12j) reduce to the first terms only and the 
two equations for F± decouple from each other. This is not the case in the presence of a 
nonlinear medium. 

The evolution takes place in an empty space, without real charges, so F-j-{r,t) have to 
satisfy the conditions 

V-F±(r,t)=0. (13) 

By applying gradient to both sides of ()12|1 it can easily be seen that V ■ F±{r, t) are constant 
in time and it is sufficient to impose the conditions ()13|) at time t = 0. 
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III. EVOLUTION OF EXEMPLARY VORTICES 

In the present section we would like to show how quantum effects connected with pairs 
creation, influence the evolution of vortices in the electromagnetic field. From amonsr the 

n 

configurations of vortex lines considered in j23| we have chosen two, for which the comparison 
can most easily be done and the effects are clearly visible. They are the situations presented 
in Figures 1 and 2 of [27]: the motion of the vortex ring and the creation and further 
evolution of initially linear vortex-antivortex configuration, i.e. two vortices of opposite 
whirl. 

Vortex lines in quantum mechanics are usually defined by the behaviour of the wave 
function of the system. In hydrodynamics vortices appear in the regions of space where 
V X v{r,t) ^ 0, where v{r,t) is the local fluid velocity. In QM, in its hydrodynamic 
formulation the role of the fluid is played by the distrubution of probability. The 
velocity field, being proportional to the gradient of the phase of the wave function, can have 
nonvanishing curl only where this phase is singular. This in turn means the vanishing of 
the wave function, i.e. the simultaneous vanishing of its real and imaginary parts. In that 
way we are led to the conclusion that, in general, vortices have the character of the curves 
(evolving in time) costituting the intersection of two surfaces defined by the requirement 
^(r,t) = 0. 

As it was proposed in j27| one can introduce in electrodynamics, in place of ip, a similar 
object, the vanishing of which may serve as the definition for the vortex lines. This object is 
F{r, ty, where F{r, t) = ^ {D{r, t) + iB{r, t)). As argued Q the quantity F is worthy of 
being called a "photon wave function" . In the case considered in the present work, photons 
move in the polarizable medium, and what is more, the nonlinear one. Already in the linear 
(but inhomogeneous) medium one is forced to define the wave function as an extention of 
F through the introduction of upper and lower components [2"^ defined by in 

nr,t)^(^*^--A, (14) 
\F^ir,t) J 

This allows one to give the set of coupled equations the form of one, linear, Schrodinger-type 
equation for the wave function J-'{r,t). In the quantum case the linearity is inevitably lost, 
but the definition of vortex lines, by the requirements S{r, t) = and V{r, t) = 0, seems 
to be universal (following 27[ this kind of singular lines has recently been called "Riemann- 



Silberstein" vortices j35|, |36[). Therefore, in the full analogy with [27[, we choose as a basic 
object the quantity 

Fl = ^{D^ - B^)+tD- B . (15) 
The condition F^{r,t)'^ = is naturally equivalent to the choice F^{r,t)'^ = 0. 



A. Vortex ring 



The first configuration considered in 27[ is defined by 

fi'')(r,t) = {y + it,z-a + i{a + t),x + it) . (16) 

This "wave function" satisfies the Maxwell equations and describes the evolution of a 
single vortex in the form of a swinging ring with varying radius. In order to see in an easy 
way how quantum (nonlinear) terms in (fT^ influence this evolution, we will choose the 
solution F^{r, t) of fT^ which is identical to (fTBj) at t = 0. This solution (up to a^) has the 
form 

(r, t) = (r, t)+t'- a{r) + f ■ /3(r) + t ■ 7(r) , (17) 

where vector funtions Q:(r), /3(r) and 7(r) are given by 

, , 128za^ , , , , 

m = ^^-(.-a)--,--a,-(a-.) + -x--a,--x + -,J , (18b) 

l{r) = ^ Qa(. - a) - ^(lla^ - 12a. + 6z^), -|a:^ . (18c) 

In Figure ^ we show the evolution of the vortex line constituting the intersection of two 
surfaces ^Fi^\r,tY = and QFP{r,tY = with Fi^\r,t) defined by (HH) and (jTHI). 
For simplicity both parameters a and m are set equal to unity on this, as well as on 



the following plots. The evo^ 
identical character to that of 



ution extends in time from t = — 1.8 to t = 1.5 and exhibits 
: the swinging vortex ring, preserving its circular character. 



decreases to certain minimal value of radius, and then starts to increase. Quantum effects 
do not manifest themselves in this domain of space and time. In the classical case, hovever, 
the expansion of a ring will last forever, which can easily be seen from the two equations 
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FIG. 1: The evolution of the vortex ring starting from t = —1.8 to t = 1.5. The scale on the axes 
is such that the frame covers the region —4 < x,y, z < 4. 

given in j27| 

x'^ + y'^ + {z- -a^ - 2at - = , (19a) 
2az + 2t{x + y + z~a)-2a^ = . (19b) 

The former represents the sphere of a fixed center in the point (0, 0, a) and of constantly 
increasing radius (for positive t). The latter, rewritten in the form x + y + {1 + a/t)z = 
a{l + a/t), tends to the motionless plane x + y + z = a passing through the center of the 
sphere. Their intersection will surely be the expanding ring. As we see in Figure |21 this 
ceases to be true in the quantum case. 

Due to the nonlinearity introduced by quantum effects two new phenomena appear. 
Firstly the vortex ring starts to deviate, for intermediate times, from its regular, circu- 
lar character. Secondly it is no longer constantly expanding. On the contrary, after reaching 
certain maximal extention it starts to decrease down to its complete disappearance, if we 
drew also frames for larger times. 

If we traced the vortex evolution even further in time (certainly far beyond the appli- 
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FIG. 2: The evolution of the vortex ring for larger times: from t = 1.5 to t = 99.3. The first frame 
is identical to the last one of Figure^ but now the scale of the axes is modified to —100 < x,y, z < 
100. 

cability of the perturbative methods) we would observe the comphcated system of vortices 
approaching from "infinity" . 

B. Vortex- ant i vortex 

The second situation corresponds to the case (b) of j^] 

/('')(r, t) = {y + t, a - i{z + a - t), X + it) . (20) 

The above function was shown to describe the configuration of two vortices which initially 
are antiparallel straight lines (we call them vortex and antivortex). They are born at t = a 
and then they separate and deform. The solution of (jl2aj) which is identical to /*^^^ at t = 
has (again up to a^) the form similar to (fTTj) 

Fi'\r,t) = f''\r,t)+t^-a{r)+t^-/3{r)+t-j{r) , (21) 
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but now with 



= -TT-I 2,— ,— , (22a) 



15m4 V ' 9 ' 9 

= (2(2; - y + a) + -a, -(x - a) + 2i{z + a), -2x + 2iy\ , (22b) 

l{r) = f ya^ + Aaz + 2^^ + + a,), 2ix\ 2iy^\ . (22c) 

The evolution of this vortex configuration is presented in Figure El Again we have put a = 1 
and m = 1. 




FIG. 3: The evolution of the system of two "antiparallel" vortices for time between t = \ and 
t = 1.55. The units on the axes are such that each frame represents the cube —4 < x, y, z < 4. 



We see in general the same motion as that found in |27l| except one difference visible in 
the first frame. In [27[ the two straight, antiparallel vortices spring up at t = a (this time 
corresponds to the first frame) as exactly overlapping. No vortices exist for t between —a 
and a. In the case of Figure El the vortices in the first frame a slightly shifted and of different 
slope. This is a result of the influence of the nonlinear (quantum) terms in (^21)j3|- 
recall that the vectors F^^\r,t) and f'^^\r,t) are synchronized for t = and not for t = a. 
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This small shift and deformation are then consequences of the quantum correction to the 
evolution for < t < a. 

In the system of vortices is born at t = a, but in our case they do not overlap and 
consequently must have appeared earlier. It is therefore interesting to take a step back in 
time and see how these vortices arise in the quantum case. This is shown in Figure IH 




FIG. 4: The appearance of the system of vortices of Figure 01 The frames correspond to times just 
before t = a. Now the scale on the y axis is changed to make the spUtting of vortices easily visible: 
-4 < X, z < 4 and -1.5 <y < 0.5. 

Figure m brings to light the essential change: the two independent vortices in the classical 
case, or rather vortex and antivortex, become the two fractions of the same, tightly bent, 
vortex line, when quantum corrections are taken into account. Their sudden creation turns 
out now to be a motion during which this single vortex line simply enters into the observation 
region and is being deformed. One might expect this kind of effects — that could be called 
topological effects — together with the smoothing of the evolution, to be the most typical 
ones introduced by the nonlinearity of the quantum equations. To make the effect more 
visible we present it again on Figure El now seen from another viewpoint. 




FIG. 5: The appearance of the system of vortices of Figure IHl seen now from the viewpoint other 
than that of Figure Eland with y-axis rescaled even more. 

We would like also to emphasize that the above phenomena take place for electromagnetic 
fields weak enough to remain in full agreement with the use of perturbation theory. 
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Yet another difference, however not visible on Figures 0] and is the shght deviation of 
the system of vortex hnes from their planar character. In the clasical case the vortex lines 
arose as the intersection of a plane with a certain surface, and therefore all vortices have to 
lie forever in one plane. This is no longer true in the quantum case. 

IV. SUMMARY 

In the present paper we analyzed the influence of the nonlinear, quantum terms in the 
Maxwell equations on the evolution of vortex lines. By making the comparison with the re- 
sults obtained earlier in the classical case {2^ we found that this evolution may be changed in 
a visible and essential way. In the first considered configuration of the constantly expanding 
vortex ring, our calculations show that quantum corrections may lead to deformation and to 
disappearance of this ring. In the second case of two linear and antiparallel vortices of the 
infinite size, which are suddenly created, we show how the process of this "creation" looks 
like, and that the two independent vortices (in the classical case) turn out to be just different 
fractions of the same vortex curve. This kind of topological changes might be expected as a 
result of nonlinearity introduced by vacuum polarization. 

The present analysis has certain limitations which come both from its perturbative char- 
acter and from the "low frequency" approximation which allowed one to derive the E-H 
Lagrangian. It can, however, serve as a qualitative picture of what type of phenomena 
may be introduced by the quantum effects. One is still very far from constructing the 
nonperturbative solutions of Quantum Electrodynamics, which would be deprived of the 
above limitations, and therefore it might be also interesting to consider the evolution of 
nodial lines in certain exact nonlinear theory as Born-Infeld electrodynamics. However, in 
this case, one cannot expect to find the polynomial solutions as given by ()17p2Hl and only 
numerical calculations come into play. This situation is similar to that in the nonlinear 
quantum mechanics. 

At the end we would like to note that although the observed deformation and evolution 
of vortices have their roots in the quantum nature of the vacuum, similar structures may 
also appear in classical and linear fields by the appropriate perturbation of the vortex con- 
figurations. Both the deviation of a vortex ring from the planar character as well as the 
occurrence of a 'hairpin'-shaped vortex, similar to that of Figure El are known in optical 
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